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Background: Phase-equivalent transformations (PETs) are well-known in quantum scattering and inverse scat- 
tering theory. PETs do not affect scattering phase shifts and bound state energies of two-body system but are 
conventionally supposed to modify two-body bound state observables such as the rms radius and electromagnetic 
moments. 

Purpose: In order to preserve all bound state observables, we propose a new particular case of PETs, a deuteron- 
equivalent transformation (DET-PET) , which leaves unchanged not only scattering phase shifts and bound state 
(deuteron) binding energy but also the bound state wave function. 

Methods: The construction of DET-PET is discussed; equations defining the simplest DET-PETs are derived. 
We apply these simplest DET-PETs to the JISP16 AW interaction and use the transformed NN interactions 
in calculations of 3 H and 4 He binding energies in the No-core Full Configuration (NCFC) approach based on 
extrapolations of the No-core Shell Model (NCSM) basis space results to the infinite basis space. 

Results: We demonstrate the DET-PET modification of the np scattering wave functions and study the DET- 
PET manifestation in the binding energies of 3 H and 4 He nuclei and their correlation (Tjon line). 

Conclusions: It is shown that some DET-PETs generate modifications of the central component while the others 
modify the tensor component of the NN interaction. DET-PETs are able to modify significantly the np scattering 
wave functions and hence the off-shell properties of the NN interaction. DET-PETs give rise to significant changes 
in the binding energies of 3 H (in the range of approximately 1.5 MeV) and 4 He (in the range of more than 9 
MeV) and are able to modify the correlation patterns of binding energies of these nuclei. 

PACS numbers: 03.65.Nk, 21.30.-x, 21.10.Dr, 21.45.Bc, 21.45.-v, 21.60.Dc 



I. INTRODUCTION 

^ ■ Phase-equivalent transformations (PETs) of two-body 
interactions are well-known in quantum scattering theory 
[l[ . PETs play an important role in the inverse scattering 
theory giving rise to ambiguities in the derived potentials. 
Currently, there is intensive research on supersymmetric 
transformations (see reviews @, H|) which are a particu- 
lar case of PETs [3| associated with removing or adding 
bound states to the system. 

More traditional PETs which do not change the on- 
shell properties of the two-body interaction, i. e. two- 
body scattering phase shifts and the energies of the two- 
body discrete spectrum states, but modify the interaction 
off-shell, were used to study manifestations of off-shell 
properties of two-nucleon interactions in many-nucleon 
systems. For example, a correlation between the nuclear 
matter binding energy and its equilibrium density (the 
so-called Coester line) was studied with phase-equivalent 
NN interactions in Ref . 5] . PETs were used to modify 



the nucleon-cluster interaction in order to obtain a cor- 
rect description of the nuclear binding energies in cluster 
model studies of Ref. @, 0] ■ Various versions of the real- 
istic JISP NN interaction (JISP6 || and JISP16 HIES]) 
were obtained by means of PETs applied to the initial 
ISTP NN interaction [ll[ obtained in the J-matrix in- 
verse scattering approach with the aim of improving the 
description of binding ener gies of many-nucleon systems. 
The interaction JISP 16 IS, 1C I provides an accurate de- 
scription of light nuclei [{J [l3-l21| and was used to pre- 
dict the binding energy and spectrum of the exotic 14 F 
nucleus |22| which were confirmed later in the first ex- 
perimental observation of this isotope (23j . 

We propose here a new type of PET, a deuteron- 
equivalent transformation (DET-PET) . Contrary to con- 
ventional PETs resulting in the modification of bound 
state and scattering wave functions 0, 0, S EH i DET- 
PET guarantees that the transformed interaction gener- 
ates not only the same scattering phase shifts and two- 
body binding energy (or, more generally, bound state 
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energies) but also the same bound state (deuteron) wave 
function as the initial untransformed interaction. The 
same method easily generalizes to preserve a set of bound 
state wave functions. DET-PET has the advantage of 
preserving the deuteron ground state observables. On 
the other hand, DET-PET, as well as any PET, modi- 
fies a two-body interaction off-shell, and hence manifests 
itself in many-body systems. 

One may naturally inquire whether PETs may lead to 
a better understanding of the appropriate off-shell be- 
havior for the NN interaction. We note that what is ap- 
propriate depends on the adopted theoretical framework 
for the NN interaction. Since the interaction is not an 
observable, all approaches (meson exchange, EFT, lattice 
gauge, inverse scattering, ...) build in model assumptions 
(e. g., form factors, regulators, cutoffs, ...). Given those 
model assumptions, there are additional unexplored off- 
shell freedoms and we will show below how to explore 
those freedoms with constraints tied to NN bound state 
observables. 

After introducing the formulas defining DET-PET, we 
apply DET-PET to the JISP16 NN interaction and illus- 
trate various versions of DET-PET by respective modifi- 
cations of scattering wave functions at a few values of the 
DET-PET continuous parameter. A DET-PET manifes- 
tation in many-body systems is illustrated by the study 
of binding energies of 3 H and He binding energies and 
their correlation (the so-called Tjon line |24|). 

It is known [H| that when any PET, DET-PET in par- 
ticular, is applied to NN interaction, the binding energy 
of a three-body (or heavier) system can be restored by 
additional three-nucleon NNN (or higher-order) inter- 
action^). Our initial 3 H applications reveal the residual 
role of the NNN interaction for the ground state energy 
and how that role changes with the DET-PET selected. 
Similarly, our initial 4 He applications reveal the residual 
roles for the combined NNN and NNNN interactions 
on the ground state energy. Given the numerical chal- 
lenges of treating NNN, NNNN, etc., interactions in 
many-body applications, it is natural to try to minimize 
their effects. In this context, DET-PETs are a poten- 
tially useful tool in future searches for an NN interaction 
consistent with many-body data. 



II. DET-PET TRANSFORMATION 

Two types of PETs are known in scattering theory: 
local PETs [l| that transform a local potential into an- 
other local potential and nonlocal PETs [HJ which gen- 
erate nonlocal potential terms. The local PETs always 
result in some modification of bound state wave functions 
[l|, IH . Therefore we focus the discussion here on nonlocal 
PETs. 

The Schrodinger equation 

H = E |* B > (1) 
describes a relative motion in two-body quantum system. 



The state \^e) can be expanded in infinite series of J;? 2 
states \a n ), 



\Ve) = ^2c n (E) \a n ). 



(2) 



n=0 



The states \a n ) are supposed to form a complete or- 
thonormalizcd basis, 



(3) 



Using expansion ^ we obtain an infinite set of algebraic 
equations defining the expansion coefficients c n (E), 



n'=0 



{H nn ,-8 nn ,E) c n ,{E)={), 



(4) 



where H nn i — (a n \H\a n i) are the Hamiltonian matrix 
elements. 

A Hamiltonian H phase-equivalent to H can be de- 
fined through its matrix [H] in the basis {\a n )}. This 
matrix [H] can be obtained from [H], the matrix of the 
Hamiltonian H in the basis {|a„)}, by means of a unitary 
transformation. 



[H] = [U}[H}[tf]. 



(•5) 



The infinite unitary matrix [U] is supposed to be of the 
form 



[U] = [U°] © [I] = 



[U°] 
[2] 



(0) 



where [I] is an infinite unit matrix and [U°] , a non-trivial 
submatrix of [U] , is a finite matrix mixing only a few se- 
lected basis functions. It is clear that Hamiltonians H 
and H have identical eigenvalue spectra. Their eigen- 
states \*S?e) and \^e), differ by a linear combination of a 
finite number of «£? 2 basis states. Any superposition of a 
finite number of J£ 2 functions must decrease at large dis- 
tances. Therefore at positive energy E associated with 
scattering, the oscillating asymptotics of wave functions 
(f\^>E) and (r\^E) at large distances are the same. In 
other words, the scattering phase shifts defined through 
asymptotic behavior of functions (t^^e) and (f\^E) are 
also the same, i. e. the Hamiltonians H and H are phase- 
equivalent. 

The unitary operator U° can be written as 



U° = 



(7) 



The transformation (|5])-([7|) leaves the bound state \d) 
unchanged, i. e. becomes a DET-PET, when each of the 
<5? 2 vectors jaj entering the non-trivial submatrix [U°] 
of the infinite unitary matrix [U] through Eq. |(7J| , is or- 
thogonal to \d), 



(ai\d) =0, i<N. 



(8) 
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At this stage, we assert that we have obtained our 
DET-PET defined through the unitary transformation 
©-(GO) w ith vectors } fitting the conditions ((3|) and 
(J5J). In order to obtain a nonlocal interaction V deuteron- 
equivalent and phase-equivalent to the initial interaction 
V, we add to V the two-body relative kinetic energy op- 
erator T to obtain the Hamiltonian H, 



H = T+V, 



(9) 



calculate its matrix [H] in the basis {|a n )}, obtain the 
matrix [H] by means of DET-PET unitary transforma- 
tion, and obtain the matrix 



[V] = [H] - [T] 



(10) 



Here [T] is the infinite kinetic energy matrix in the basis 
{|a n )}. The interaction V is defined through its matrix 
[V] in the basis {|a n }}. 

The simplest DET-PET is obtained with arbitrary uni- 
tary matrix [U ] of the rank 2. In this case, [U°] is 
associated either with a rotation by the angle j3 when 
dct U° — +1 or with a rotation by the angle /3 combined 
with reflection when det U° = — 1 . We also need to define 
the J? 2 vectors |oi) and \a 2 ) in Eq. 0. 

We define here the vectors |ai) and \a 2 ) as linear com- 
binations of oscillator states | 



i«*> = E 



(ii) 



'<N' 



which fit the orthonormality condition ([3]). We expand 
the deuteron eigenstate \d) in an infinite series of oscilla- 
tor states, 



(12) 



i=0 



where, generally, all the coefficients di are non-zero, 

dj ¥= 0. (13) 

Since the vectors |oi) and \a 2 ) should fit Eq. ©, the 
expansion (fTTj) of each of them involves at least two dif- 
ferent basis states \(fi) due to Eq. P^|) - (fT3^) . In this 
simplest case we have 



\ai) = a? \ip n ) + a™ \ip m ), 
\a 2 ) = a 2 \ip k ) +4 \ipi). 



(14a) 
(14b) 



The normalization of these vectors requires 

«) 2 + (a?) 2 = 1, (15a) 
(a 2 fe ) 2 + (4) 2 = l, (15b) 

while the orthogonality of the vectors \ai) and \a 2 ), 

(a 2 \a 1 )=0, (16) 



is guaranteed when these vectors are constructed from 
different basis states, i. e. all the basis states \<p n ), \f m ), 
\<Pk)t \<Pi) entering Eqs. lfT4")) are different. Using expan- 
sions (fTTj) and (fl"2"|) we obtain 



a? d n +a 7 [ l d m = Q, 
a 2 dk + a l 2 di = 0. 

The solutions of Eqs. {TSJ, HZ]) are 

a „ _ d m 
m d n 



(17a) 
(17b) 



(18a) 
(18b) 
(18c) 
(18d) 



To define completely the simplest DET-PET discussed 
above we need to fix the rotation angle (3, the sign of 
det U° and the set of 4 oscillator states used to build the 
states | ax) and \a 2 ). To distinguish various DET-PET 
types we use notations like 0s2sls2d ± . In this example, 
the state \ai) is a linear combination of the oscillator 
states 0s and 2s, the vector \a 2 ) is a linear combination 
of the oscillator states Is and 2d, and the index ± corre- 
sponds to the sign of det U° = ±1. 



III. DET-PET PROPERTIES AND 
MANIFESTATION IN FEW-NUCLEON SYSTEMS 

In this section, we study modifications of the JISP16 
NN interaction @ induced by various DET-PETs. The 
modifications of a nonlocal interaction can be illustrated 
by modifications of its wave functions. The deuteron 
wave function is unaffected by DET-PET. Therefore we 
present below the DET-PET induced transformation of 
the JISP16 np scattering wave function in the sd coupled 
partial wave. 

It is interesting to explore a DET-PET which acts only 
in a single channel, say, in the s channel, and compare it 
with DET-PETs mixing components of the s and d chan- 
nels in different ways. Therefore vectors |ai) and \a 2 ) [see 
Eqs. (fT4)) ] were constructed as various superpositions of 
two low-lying oscillator states of the np relative motion 
0s, Is, 2s, 3s, Od and Id with Ml = 40 MeV. For each 
type of the DET-PET we investigate the transforma- 
tions associated with both pure rotation and a rotation- 
reflection combination. 

Plots of the JISP16 np scattering wave functions in 
the sd coupled partial wave at laboratory energy E\ & \, = 
10 MeV and plots for their 0s2sls3s ± , OslsOdlc^ and 
lsOdOsl^ DET-PET partners are given in Figs. [2 H 
and |3] respectively. We use the i^T-matrix formalism (see 
Ref. [ll| for details and nomenclature adopted here). The 
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FIG. 1. (Color online) Large (a) and small (b) components of the np scattering wave function at the laboratory energy 
Slab = 10 MeV in the sd coupled partial wave in the _R"-matrix formalism (see Ref. [Tl| f° r details and nomenclature) generated 
by JISP16 and NN interactions obtained from JISP16 by means of DET-PET 0s2sls3s ± . The sign of det U° is given in the 
legends in parenthesis after the value of rotation angle ft. 




FIG. 2. (Color online) Same as Fig. [T]but for DET-PET Qsls0dld ± . 



advantage of the fC-matrix formalism is that the radial 
wave functions in the scattering domain defined accord- 
ing to their standing wave asymptotics are real contrary 
to the more conventional S'-matrix formalism with com- 
plex radial wave functions which are asymptotically a 
superposition of ingoing and outgoing spherical waves. 

The DET-PET 0s2sls3s ± mixes various s wave com- 
ponents of the wave function that is equivalent to modifi- 
cation of the central part of the JISP16 interaction in the 
s wave. This results in significant changes of the large s 
wave component as is seen in Fig. [TJ The modification of 
the small s wave component is less pronounced. The d 
wave components, as expected, are nearly unaffected by 
0s2sls3s ± . 

The DET-PET OslsOdld* explicitly mixes s and d 



waves; the DET-PET lsOdOsld also mixes s and d waves 
but in a different manner. This corresponds to an essen- 
tial modification of the tensor component of the JISP16 
NN interaction. As a result, we observe an essential 
modification of small scattering wave function compo- 
nents which are generated by the tensor NN interaction 
as is seen in Figs. [5] and [3] Modifications of the large 
wave function components are much less pronounced. 

We see that DET-PET generates essential modifica- 
tions of scattering wave functions without any change of 
scattering phase shifts and scattering observables (cross 
sections, polarization observables, etc.). It is worth not- 
ing here that the deuteron wave function and deuteron 
observables (rms radius, quadrupole moment, etc.) are 
unaffected by DET-PET due to the nature of this trans- 
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formation. The alteration of scattering wave functions is 
an indicator for the variation of the off-shell properties 
of the interaction arising from DET-PET. The modifica- 
tion of the the NN interaction off-shell should manifest 
itself in properties of many-nucleon systems. Therefore, 
we investigate the DET-PET-induced changes of the 3 H 
and 4 He binding energies. 

We calculate 3 H and 4 Hc in the ab initio No-core Full 
Configuration (NCFC) [13 approach. Within the NCFC 
approach, we start with the No-core Shell Model (2(| \2l\ 
calculations using the code MFDn 28-3jJ with a few 
values of the oscillator frequency Ml and in a few basis 
spaces characterized by the maximum oscillator quanta 
iV max allowed in the many-body basis above the minimal 
configuration. Next, we extrapolate the sequence of finite 
basis space results to the infinite basis space limit. This 
makes it possible to obtain basis space independent re- 
sults for binding energies and to evaluate their numerical 
uncertainties. NCFC suggests two extrapolation meth- 
ods: a global extrapolation based on the calculations 
in four successive basis spaces and five Ml values in a 
10 MeV interval (extrapolation A), and extrapolation B 
based on the calculations at various fixed Ml values in 
three successive basis spaces and defining the most re- 
liable Ml value for the extrapolation. We present here 
only the extrapolation A results based on the NCSM cal- 
culations with basis spaces up through N max = 16. The 
extrapolations A and B usually provide consistent results 
[17) . and we checked this consistency for our results in a 
number of cases. The evaluated uncertainties of results 
for binding energies presented here are less then 10 keV 
in most cases; in a few cases, we performed the NCSM 
calculations up to -ZV max = 18 to obtain the binding en- 
ergies with uncertainty of about 10 keV. 

The binding energies of 3 H nucleus E t and of He 
nucleus E a were calculated with JISP16 interaction 
modified by DET-PETs 0s2sls3s ± , 0sls0dld ± and 
Is0d0sld ± varying angle /3 from 0° through 360° in steps 



of 60° . We observe variations of E t and E a due to DET- 
PETs. In some cases, when the 3 H and 4 He binding 
energies were close to their maximal or minimal values 
for a given DET-PET type, we decreased the step of /3 
to investigate the behavior of E t and E a around their 
extremal values in more detail. 

The ranges of 3 H and 4 He binding energy variations 
for each DET-PET type are shown in Table [I] We see 
that DET-PETs can cause essential modification of both 
3 H and 4 He binding energies. For example, in the case of 
the 4 He nucleus, E a can be varied by DET-PETs on the 
interval from 21.25 through 30.41 MeV, i. e., the DET- 
PET NN interaction can change E a by more than 7 MeV 
from its original value provided by the original JISP16 
interaction. In the case of 3 H, the range of the DET-PET 
binding energy variation is 7.21 < E t < 8.67 MeV, i. e., 
the binding energy can be shifted by more than 1 MeV 
from its original JISP16 value. 

We study also a correlation of the 3 H and 4 He bind- 
ing energies, the so-called Tjon line (24|. The Tjon line 



TABLE I. Ranges of 3 H and 4 He binding energy variations 
(in MeV) caused by various types of DET-PET in compari- 
son with the binding energies obtained with JISP16 and their 
experimental values. 



'He 



'He 



0s2sls3s+ 
7.21—8.37 21.25—28.49 

0sls0dld+ 
7.67—8.41 23.50—28.83 

Is0d0sld+ 
7.98—8.64 25.79—30.36 
JISP16 
8.369(1) 28.299(1) 



0s2sls3s 
7.25—8.35 21.46—28.59 

OslsOdld' 
7.68—8.39 23.46—28.91 

lsOdOsld" 
8.05—8.67 26.18—30.41 
Experiment 
8.482 28.296 
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FIG. 4. (Color online) Tjon line obtained with DET-PET 
Is0d0sld ± in comparison with results obtained with various 
NN and NN + NNN interaction models from Refs. [H, SJ, 



is usually studied using results obtained with different 
NN interactions and different combinations of iV/V and 
NNN interactions (see, e. g., Ref. [3211 ). We note here 
also an investigation of Jurgenson et al [33| where the 3 H 
and 4 He binding energy correlation was studied with 7Y7Y 
interactions SRG evolved to various values of momen- 
tum parameter A. An interesting observation mentioned 
by various authors (see, e. g., Ref. [HlHij]) is that these 
results obtained with different interaction models form 
nearly a straight line on the plot E a vs E t . Here we study 
the E t -E a correlation using families of NN potentials 
generated by various DET-PET types from the JISP16 
interaction, i. e., all NN interactions provide not only 
algebraically identical NN phase shifts but also identical 
deuteron wave functions that should give rise to specific 
np correlations in three- and four-nucleon systems. 

We begin the discussion of the Tjon lines from the 
results obtained with the Is0d0sld ± DET-PET pre- 
sented in Fig. [4] where we show also the results from 
Refs. [iH, [34], l3a | obtained with various modern NN 
and NN + NNN interaction models. It is seen that the 
DET-PET E t -E a correlation generally follows the trend 
suggested by other interactions: our results are concen- 
trated close to the Tjon line connecting the points ex- 
tracted from Refs. [HI, [34], HU and extend it to larger 
3 H and 4 He bindings. We recall here that the lsOdOsld^ 
DET-PET is associated with modification of the tensor 
component of NN interaction. 

Another DET-PET modifying the tensor component of 
the JISP16 NN interaction, is the DET-PET 0sls0dld ± . 
This DET-PET results in a very different range of E t 
and E a variations (see Table H}. The DET-PET 3 H and 
4 He binding energies are also correlated along a nearly 
straight line (see Fig. [5]). However this line has a slope 
very different from the slope of the Tjon line obtained 
with other interaction models. Around the maximal 3 H 



" *-»0sls0dld 
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FIG. 5. (Color online) Same as Fig. H but for DET-PET 



•-•0s2sls3s 




E t (MeV) 

FIG. 6. (Color online) Same as Fig. H but for DET-PET 
0s2sls3s ± . 



and 4 He binding energies accessible by this DET-PET, it 
suggests correlations consistent with those derived using 
modern NN + NNN interaction models. However, for 
smaller binding energies, this DET-PET suggests much 
less bound 4 He at the same 3 H bindings as provided by 
modern purely two-nucleon interactions. 

The DET-PET C^slsSs* modifies the central s-wave 
component of the NN interaction. It results in the 3 H 
and 4 He binding energy correlation shown in Fig. [6] We 
see that in this case the results do not concentrate as 
tightly around some straight line. That is, they are more 
spread out on the E t -E a plane. The DET-PET Tjon 
lines transform into closed-loop curves surrounding elon- 
gated areas. In the case of the DET-PET 0s2sls3s~, the 
Tjon curve surrounds many points obtained with various 
NN interactions. The DET-PET 0s2sls3s+ generates 
the Tjon curve shifted down from the Tjon line suggested 
by other interactions. Both 0s2sls3s + and 0s2sls3s~ 
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DET-PETs essentially extend the range of the 3 H and 
4 He binding energy variations to smaller bindings. 



IV. CONCLUSIONS 

We have introduced a new type of phase-equivalent 
transformations, DET-PET, preserving the deuteron 
wave function. The suggested theory of DET-PET can 
be easily reformulated to preserve scattering wave func- 
tions at a given energy instead of the bound state wave 
function. We investigated transformations of the JISP16 
NN interaction induced by DET-PETs mixing oscilla- 
tor components in various combinations. One of these 
DET-PETs generates modifications of the central com- 
ponent of the NN interaction, the others modify the ten- 
sor NN interaction component. We demonstrated that 
DET-PETs are able to modify significantly the np scat- 
tering wave functions and hence the off-shell properties of 
the NN interaction while the on-shell interaction prop- 
erties are preserved. 

DET-PETs impact the binding energies of many- 
nucleon systems. We found that the 3 H and 4 He binding 
energies can be significantly changed by DET-PETs. The 
investigated DET-PETs modifying tensor NN interac- 
tion, correlate the 3 H and 4 He bindings along some lines 
that may differ in slope from the Tjon line obtained with 
modern NN and NNN interactions. The DET-PET 
0s2sls3s ± modifying the central s wave NN interaction, 
weakens the E t -E a correlation spreading the results on 
the E t -E a plane. 



It would be interesting to study DET-PET manifesta- 
tions in binding energies and other observables of heavier 
nuclei. We speculate that DET-PET can be helpful in the 
further development of JlSP-like NN interactions. 

DET-PETs can be also used to design an interesting 
approach to effective interactions. In particular, DET- 
PETs can be applied to a modern NN realistic interac- 
tion to reduce their high momentum components (and 
hence to improve the convergence of ab initio calcula- 
tions). At the same time, we are preserving the deuteron 
wave function and NN correlations in other partial waves 
at some energy by using the extension of DET-PET to 
preserve the scattering wave function. Such an effective 
interaction can be very interesting for many-body nuclear 
applications. 

It is also possible to extend the DET-PET concept to 
NNN interactions. The corresponding set of transfor- 
mations would then involve changing the interior part of 
the NNN wave function in such a manner as to preserve 
the NNN ground state wave function and to preserve 
the asymptotic NNN wave function. This could provide 
a useful tool to explore the off-shell freedoms available in 
NNN interactions without requiring repeated fits to the 
NNN bound state properties. 

We are thankful to Peter Sauer for valuable dis- 
cussions. This work was supported by the US DOE 
Grants DE-FC02-09ER41582 and DE-FG02-87ER40371 
and Ministry of Education and Science of the Russian 
Federation under the Contract P521. 
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